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Abstract 

We make use of product integrals to provide an unambiguous mathematical 
representation of Wilson line and Wilson loop operators. Then, drawing upon 
various properties of product integrals, we discuss such properties of these op- 
erators as approximating them with partial sums, their convergence, and their 
behavior under gauge transformations. We also obtain a surface product inte- 
gral representation for the Wilson loop operator. The result can be interpreted 
as the non-abelian version of Stokes theorem. 

1 Introduction 

The notion of Wilson loop U provides a systematic method of obtaining gauge 
invariant observables in gauge theories. Its applications range over such diverse fields 
as phenomenology and lattice gauge theories on the practical side and topological 
gauge theories || and string theory |3J on the purely theoretical side. The impor- 
tance of the Wilson line as a parallel transport operator in the gauge independent 
formulation of gauge theories has been emphasized by Mandelstam ||, and further 
developed by Wu and Yang ||. More recently, in the context of the AdS/CFT cor- 
respondence [0], an interesting connection between Wilson loops in supersymmetric 
gauge theories and membranes in supergravity theories has been suggested M. In 
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view of all these developments, it is imperative that Wilson lines and Wilson loops 
be described within a well defined mathematical framework. The main purpose of 
this work is to provide such a representation by means of product integrals. This will 
permit us to give, among other things, two unambiguous proofs of the non-abelian 
version of Stokes theorem. 

The product integral formalism has been used extensively in the theory of differ- 
ential equations and of matrix valued functions |J. In the latter context, it has a 
built-in feature for keeping track of the order of the matrix valued functions involved. 
As a result, product integrals are ideally suited for the description of path ordered 
quantities such as Wilson lines and Wilson loops. In fact, they make precise what 
one means by these concepts as well as what one means by the notion of path or- 
dering in general. Moreover, since the theory of product integrals is well developed 
independently of particular applications, we can be confident that the properties of 
Wilson lines and Wilson loops which we establish using this method are correct and 
unambiguous. 

Among the important advantages of the product integral representation of Wilson 
lines, one is the manner in which it deals with convergence issues. In the physics 
literature, the exponential of an operator such as a Wilson line is defined formally 
in terms of its power series expansion. In such a representation, it will be difficult, 
without any further elaboration, to establish whether the series converges and if so 
how. In contrast to this, it is a straight forward matter to establish the criteria for the 
convergence of the Wilson line in its product integral representation. This is because 
in such a framework the Banach space structure of the corresponding matrix valued 
functions is already built into the formalism. 

Another important advantage of the product integral formulation of the Wilson 
loop is that, at least for orientable surfaces, it permits a 2-surface representation for 
it. Based on the central role of Stokes theorem in physics and in mathematics, it 
is not surprising that the non-abelian version of this theorem has attracted a good 
deal of attention in the physics literature [|T0f] - fl20| . The central features of the earlier 
attempts[T^]- [|T^] have been reviewed and improved upon in a recent work |]I7[| . Other 
recent works on non-abelian Stokes theorem |18|, [19], focus on specific problems 
such as confinement ]nj, zig-zag symmetry |2(J suggested Polyakov [PTJ] , etc. With 



one exception JT7|], the authors of these works seem to have been unaware of a 1927 
work in the mathematical literature by Schlessinger [2^] which bear strongly on the 
content of this theorem. Modern non-abelian gauge theories did not exist at the 
time, and Schlesinger's work dealt with integrals of matrix valued functions and their 
ordering problems. Its relevance to Wilson lines and Wilson loops is tied to the fact 
that in non-abelian gauge theories, the connection and curvature are matrix valued 
functions. As a result, Schlesinger's work amounts to establishing the non-abelian 
Stokes theorem in two (target space) dimensions. By an appropriate extension and 
reinterpretation of his results, we show that the product integral approach to the 
proof of this theorem is valid in any target space dimension. 



2 



This work is organized as follows: To make this manuscript self- contained, we 
review in Section 2 the main features of product integration || and state without 
proof a number of theorems which will be used in the proof of the non-abelian Stokes 
theorem and other properties. In Section 3, we express Wilson lines and Wilson 
loops in terms of product integrals. In Section 4, we turn to the proof of the non- 
abelian Stokes theorem for orientable surfaces. In section 5, we give a variant of this 
proof. Section 6 is devoted to convergence issues for Wilson lines and Wilson loops. 
Section 7 deals with two observables from the Wilson loop operator. In Section 8, 
we study behavior of Wilson lines and Wilson loops under gauge transformations. 
This is another instance in which the significant advantage of the product integral 
representation of these operators becomes transparent. 

2 Some properties of product integrals 

The method of product integration has a long history, and its origin can be traced 
to the works of Volterra||. The justification for its name lies in the property that 
the product integral is to the product what the ordinary (additive) integral is to the 
sum. One of the most common applications of product integrals is to the solution of 
systems of linear differential equations. To see how this comes about, let us consider 
an evolution equation of the type 



where s is a real parameter, and prime indicates differentiation. When the quantities 
Y and A are ordinary functions, and Y Q is an ordinary number, the solution is given 
by an ordinary integral. On the other hand, if these quantities are matrix valued 
functions arising from a system of, say, n linear differential equations in n unknowns, 
then the solution will be a product integral. 

To motivate a more precise formulation of product integration, we start with a 
simple example which exhibits its main features. Let us suppose that all the matrix 
valued functions that appear in the above equation are continuous on the real interval 
[a, b}. Then, given the value of Y at the point a, i.e., given Y(a), we want to find 
Y(6). One can obtain an approximate value for Y(6) using a variant of Euler's 
tangent-line method. Let P = {s , si, . . . , s n } be a partition of the interval [a, b], and 
let Ask = Sk — Sfc_i for all k — 1, . . . , n. In the interval [sq, si], we approximate A(s) 
by the constant value A(s\), solve the differential equation with initial value Y(a) 
and get the approximate solution for Y at s^: 



In the next interval [si,s 2 ], using the above approximate value as input, and 
replacing A(s) by Afa), one finds 
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Proceeding in this manner we obtain the following approximate value for Y(6): 

Y(6) « e A(s " )As " . . . e A(si)Asi Y(a) = U p (A)Y(a) : (2) 

where 

n 

H p (A) = J] e A(Sk)Ask . 

k=l 

We stress that the order of the exponentials on the right hand side of this equation 
is important since the corresponding matrices do not commute in general. 

Since A is continuous on the compact interval [a, b], it follows that A is uniformly 
continuous in that interval. This implies that for all k = 1, . . . , n the value A(sj.) will 
be close to the values of A(s) on [sfe_i, sjt]. It is thus reasonable to suppose that if 
the mesh fx(P) of the partition P (the length of the longest subinterval) is small, the 
above calculation results in a good approximation to Y(6). Then, we expect that the 
exact value of Y(6) is given by the natural limiting procedure 

Y(6) = ^mU P (A)Y(a) ee f{e A ^ ds Y(a). (3) 

Having identified the main ingredients which characterize the above construction, 
we proceed to give a precise definition of the product integral || . We begin with the 

Definition 1 Let A : [a, b] — > C nxn be a function with values in the space of complex 
nx n matrices. Let P = {sq, s±, . . . , s n } be a partition of the interval [a,b], with 
As/c = Sk — Sjfe-i for all k — 1, . . . , n. 

(i) A is called a step function iff there is a partition P such that A is constant on 

each open subinterval (sfe-i, Sk) for all k = 1, . . . , n. 

(ii) The point value approximant Ap corresponding to the function A and partition 

P is the step function taking the value A(sk) on the interval (sfc_i, Sfc] for all 
k — 1, . . . , n. 

(Hi) If A is a step function, then we define the function Ea : [a,b] — > C nxn by 
E A {x) := e M»k)(x-»k-i) . . . e A(s 2 )As 2e A(s 1 )As 1 f Qr any x e ( Sk _ 1} s k ], for all k = 

1, ... ,n, and Ea{o) := /. 

After a number of intermediate developments, one arrives at the following fundamen- 
tal theorem which can be taken as the starting point of product integration: 

Definition-Theorem 1 Given a continuous function A : [a, b] — > C nxn and a se- 
quence of step functions {A n }, which converge to A in the sense of L l ([a, b]), then 
the sequence {Ea„{x)} converges uniformly on [a,b] to a matrix called the product 
integral of A over [a,b]. 
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More explicitly, we have: 

b 

The product integral of A over [a, b] = e w*. (4) 

a 

Now we are in position to enumerate some of the basic properties of product 
integrals. The proofs of these assertions are given in reference ||. Let A : [a, b] — > 
C nxn be a continuous function, and for any x G [a, b] let 

F(x,a) :=f[e A{s)ds (5) 

a 

denote the product integral from a to x. Then, F satisfies the following integral 
equation: 

F(x,a) = 1+ f dsA(s)F(s,a), (6) 

J a 

where I = I nxn is the n x n unit matrix. The quantity F is also a solution of the 
following initial value problem: 

d 

—F(x,a)=A(x)F(x,a), F(a,a) = I. (7) 
dx 

Although product integrals can formally be defined for singular matrices, the above 
definition makes sense if they are non-singular. This is true, e.g., when the matrices 
form a group. Then the determinant of the product integral is given by the following 
theorem: 

Theorem 1 Given the continuous function A : [a, b] — > C nxn , then for every x G 
[a, b], the product integral ]la e A ^ ds is non-singular and the following formula holds: 

det (fle A ^ ds ^j = e X> AW(fc , (8) 

where "tr" stands for trace. 

When the set of matrices {A(s) : s G [a, b}} is commutative, i.e. [j4(s), A(s')] = 
Vs, s' G [a, b], it is easy to show that 

Y[ e A(s)ds = Jl A{ S )ds_ ^ 
a 

It is convenient to define the product integral Y\ b a e A( - s ^ ds also in the case when a > b. 
It will be recalled that for ordinary (additive) integrals f a A(s)ds = — J fe a A(s)ds. To 
obtain the analog of this for product integrals, we merely replace the "additive" 
property with the corresponding "multiplicative" property: 

Y[ e Ms)ds ._ (f[ e A(s)ds\ (1Q) 
a \ b ) 
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The additive property of ordinary integrals also provides a composition rule for them: 
j° A(s)ds + f h c A(s)ds = J a b A(s)ds. For product integrals, we have an analogous 
composition rule 0: 

b be 
a c a 

Another well known property of ordinary integrals is the differentiation rule with 
respect to the endpoints: ^ f/' A(s)cZs) = A(6) and £ f/ a 6 A(s)ds) = — A(a). The 
following theorem gives the corresponding rule for product integrals. 

Theorem 2 Let A : [a, 6] — > C nxn be continuous. For any x, y G [a, 6] we /iai>e: 

It is important to keep in mind the relative order of the A(x) and A(y) with respect 
to the product integral. 

The usual elementary way of computing ordinary integrals is by means of the 
fundamental theorem calculus: J a f(x)ds = F(b) — F(a), where F is a primitive 
function of / (F' — f). To obtain the corresponding theorem for product integrals, we 
start by defining the so called L-operation which is a generalization of the logarithmic 
derivative, for non-singular functions: 

Definition 2 The L-derivative of a non-singular differentiable function P : [a, b] — ► 
is given by: 

LP(x) := P'{x)P-\x). (13) 

To demonstrate the usefulness of this operation, let us consider the product integral 
p(x) = ]\ x a e A ^ ds P(a). From Theorem g, we have P'(x) = A(x) \[ x a e A ^ ds P(a) = 
A(x)P(x). Then, from the above definition, we get (LP)(x) = A(x) (the derivative 
of the primitive function is the original function). We are thus led to the analog of 
the fundamental theorem of calculus for product integrals 0: 

Theorem 3 For a non-singular and continuously differentiable (C [a, b]) function 
P : [a, b] — > C nxn , we have 

f[e {LP){s)ds = P(x)p-\a). (14) 

a 

The following elementary properties of the L-operation follow from its definition: 

{LP- l )(x) = (p-'YiP- 1 )- 1 = (-P^P'P-^P = -P-\x)P / (x), (15) 

and 

L(PQ)(x) = (P'Q + PQ')Q _1 P _1 = LP(x) + P(x)(LQ(x))P-\x). (16) 

We will rely heavily on the contents of the next three theorems in proving the non- 
abelian version of Stokes theorem. The proofs are given in reference 0. 
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Theorem 4 (Sum rule): Given continuous functions A,B : [a, b] — > C nxn , let 
P(x) = U x a e A ^ ds . Then 

Y^ e [A(s)+B(s)]ds _ p^YT e P- 1 (s)B{s)P(s)ds_ ^ 



Theorem 5 (^Similarity rulej: Given a continuous function B : [a, b] — > C nxn and 
the non-singular function P : [a, b) — > C nxn , then 



s)+P(s)B{s)p- 1 (s)}ds 



P(x) (f[e B ^ ds ) p-\a) = f[e^ 



Theorem 6 (Derivative with respect to a parameter ): Given a function A : [a, b] x 
[c,d] — > C nxn such that A(s, A) is continuous in s for each fixed X e [c,d). and is 
differentiable with respect to A. Then, the product integralP(x,y; A) = JI^ e A( - s ' ,x "> ds is 
differentiable with respect to X, and 



d f x dA 

—P(x,y;X) = J y dsP(x,s;X)—(s;X)P(s,y;X). (19) 



To put the above description in its proper perspective, we note that instead of the 
specific complex Banach space L 1 ([a, b], product integrals can be defined over more 
general Banach spaces. Consider, e.g., a set B(X) of bounded linear operators over 
a complex Banach space, and let A : [a, b] — > B(X) be an operator valued function. 
It is possible to define the product integral of A and establish the analogs of the 
properties given above in this more general context. Then, the standard topologies 
(norm, strong, and weak) on the space of bounded linear operators play an important 
role. Moreover, the notion of Lebesgue integrable functions used on L 1 ([a, b}) space 
above generalize naturally to Boschner integrable functions [23|. For details we refer 
again to 0. 



3 The representation of Wilson Line and Wilson 
Loop 

As noted in the introduction, Wilson lines and Wilson loops are intimately related 
to the structure of non-abelian gauge theories. To provide the background for using 
the product integral formalism of Section 2 to explore their physical properties, we 
begin with the statement of the problem as it arises in the physics context. Let M 
be an n-dimensional manifold representing the space-time (target space). Let A be 
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a (connection) 1-form on M. When M is a differentiable manifold, we can choose a 
local basis dx 11 , fi = 1, ...,n, and express A in terms of its components: 

A(x) = Apix) dx^. 

We take A to have values in the Lie-algebra, or a representation thereof, of a Lie 
group. Then, with Tj., k — 1, ..,m, representing the generators of the Lie group, the 
components of A can be written as 

A,(x) = Aj(x) T k . 

With these preliminaries, we can express the Wilson line of the non-abelian gauge 
theories in the form IE3 



w ab (c) = vX A . 

where V indicates path ordering, and C is a path in M. When the path C is closed, 
the corresponding Wilson line becomes a Wilson loop |23: 



W{C)=Ve$ A . (20) 

The path C in M can be described in terms of an intrinsic parameter a, so that for 
points of M which lie on the path C, x^ 1 = x^{a). One can then write 

A IJ ,(x(a))dx tl = A(a)dcr, 

where 

A{a) = A»(x(o))^l. 

da 

It is the quantity A(a), and the variations thereof, which we will identify with the 
matrix valued functions of the product integral formalism. 

Let us next consider the Wilson loop. For simplicity, we assume that M has 
trivial first homology group with integer coefficients, i.e., H\(M, Z) = 0. This insures 
that the loop may be taken to be the boundary of a two dimensional surface £ in 
M. More explicitly, we take the 2-surface to be an orientable submanifold of M. It 
will be convenient to describe the properties of the 2-surface in terms of its intrinsic 
parameters a and r or a a , a = 0, 1. So, for the points of the manifold M, which lie 
on E, we have x = x(a, r). The components of the 1-form A on E can be obtained 
by means of the vielbeins (by the standard pull-back construction): 

Thus, we get 

A a = v a A^. 

The curvature 2-form F of the connection A is given by 

F = dA + AAA = -F^ dx 11 A dx v . 



8 



The components of F on £ can again be obtained by means of the vielbeins: 

1 ab u a u b 1 W 

We note at this point that we can construct the pulled-back field strength F a ^ in 
another way, as the the field strength of the pulled-back connection A a . It is easy to 
check that these two results coincide, insuring the consistency of the construction. 

We want to express the Wilson loop operator in terms of product integrals [ p5fl . 
To achieve this, we begin with the definition of a Wilson line in terms of a product 
integral. Consider the continuous map A : [so,Si] — > C nxn where [s ,Si] is a real 
interval. Then, we define the Wilson line given above in terms of a product integral 
as follows: 

so 

Anticipating that we will identify the closed path C over which the Wilson loop is 
defined with the boundary of a 2-surface, it is convenient to work from the beginning 
with the matrix valued functions A(a,r). This means that our expression for the 
Wilson line will depend on a parameter. That is, let 

A : [<r , (Ji ] x [t , Ti] -> C nxn , (21) 

where [o"o,(7i] and [to,Ti] are real intervals on the two surface £ and hence in M. 
Then, we define a Wilson line 

P(a,a ;r) = f[ e A ^'^' ee vJ^'^' . (22) 

CO 

In this expression, V indicates path ordering with respect to a, while r is a parameter. 
To be able to describe a Wilson loop, we similarly define the Wilson line 

Q(a; r, r ) = f[ e A ^'^' = vA A ^ T ' )dr ' . (23) 

TO 

In this case, the path ordering is with respect to r, and a is a parameter. 

To prove the non-abelian version of the Stokes theorem, we want to make use 
of product integration techniques to express the Wilson loop operator as an integral 
over a two dimensional surface bounded by the corresponding loop. In terms of the 
intrinsic coordinates of such a surface, we can write the Wilson loop operator in the 
form 

W(C)=Pef A * daa , (24) 

where, as mentioned above, 

a a = (r,a); a =(0,1). (25) 
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The expression for the Wilson loop depends on the homotopy class of paths in M to 
which the closed path C belongs. We can, therefore, parameterize the path C in any 
convenient manner consistent with its homotopy class. In particular, we can break up 
the path into segments along which either a or r remains constant. The composition 
rule for product integrals given by Eq. (11) ensures that this break up of the Wilson 
loop into Wilson lines does not depend on the intermediate points on the closed path 
which are used for this purpose. So, inspired by the typical paths which are used in 
the actual computations of Wilson loops (see for example 0), we write 

W = W 4 W 3 W 2 W ll (26) 

In this expression, Wk, k = 1, ..,4, are Wilson lines such that r = const, along W\ 
and W3, and a = const, along W 2 and W4. 

To see the advantage of parameterizing the closed path in this manner, consider 
the exponent of Eq. ^ : 

A a da a = A dr + A x da. (27) 

Along each segment, one or the other of the terms on the right hand side vanishes. 
For example, along the segment [a , a], we have r = r = const.. As a result, we get 
for the Wilson lines W\ and W 2 , respectively, 

Wl = n e A ^'^' = Ve J °o A ^™W = p( t7j ao; ro ) 9 (28) 

TO 

and 

w 2 = n e M°-y)dr> = V J T T0 M«-s) d T> = Q{a _ r ro) _ (29) 

TO 

When the 2-surface S requires more than one coordinate patch to cover it, the 
connections in different coordinate patches must be related to each other in their 
overlap region by transition functions |J. Then, the description of Wilson loop in 
terms of Wilson lines given in Eq. (26) must be suitably augmented to take this 
complication into account. The product integral representation of the Wilson line 
and the composition rule for product integrals given by Eq. (11) will still make it 
possible to describe the corresponding Wilson loop as a composite product integral. 
For defmiteness, we will confine ourselves to the representation given by Eq. (26). 

It is convenient for later purposes to define two composite Wilson line operators 
U and T according to 

U(a, r) = Q(a; t, r ) P(a, a ; r) , (30) 

T(a;r) = P(a,a ]r) Q(a ;T,T ). (31) 
Using the first of these, we have 

W 2 W 1 = U(a,r). (32) 
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Similarly, we have for the two remaining Wilson lines 

W 3 = P- 1 (a, ( xo;T), (33) 

and 

W 4 = Q- l {a -T,T ). (34) 
From the Eq. (|3~TD, it follows that 

W i W 3 = T~ 1 (a,r). (35) 

Appealing again to Eq. (11) for the composition of product integrals, it is clear that 
this expression for the Wilson loop operator is independent of the choice of the point 
(er, r). In terms of the quantities T and U, the Wilson loop operator will take the 
compact form 

W = T- 1 (a;r)U(a;r). (36) 



4 Non-Abelian Stokes Theorem 

As a first step in the proof of the non-abelian Stokes theorem, we obtain the action 
of the L-derivative operator on W: 

L T W = L T [T-\a, r)Q(a; r, r )P(a, a ; r )]. (37) 

Using the definition of the L-operation given by Eq. (13), noting that P(a, a ; r ) is 
independent of r, and carrying out the L operations on the right hand side (RHS), 
we get 

L T W = L T T- 1 (a,r)+T- 1 (a,r)[L r Q(a;r,r ) + 

+Q(a; r, T )(L T P(a, a o] r ))Q-\a; r, r )]T(a, r). (38) 

Simplifying this expression by means of Eqs. (13) and (15), we end up with 

L T W = T-\a, r)[A (a, r) - L T T(a, r)]T(a, r). (39) 

Next, we prove the analog of Theorem 3, which applies to an elementary Wilson line, 
for the composite Wilson loop operator defined in Eq. (24) and made explicit in Eq. 
(26). 

Theorem 7 The Wilson loop operator defined in Eq. (26) can be expressed in the 
form 

W = fl e T ~ 1(ay)[Ao{ ° y) ~ LTT{a ' T ' )]T(cTy)dT ' . (40) 

T 



11 



To prove this theorem first we note from the definition of the L operation that the 
right hand side (RHS) of this equation can be written as 



RHS = r[e [T " 1(<7;T ' )Ao(CT;r ' )T(CT;T ' )_T " 1(<T;r,) ^ T(<T;T ' )ldT '. (41) 

TO 

Noting that —T~ 1 d T T = L T T, we can use Theorem [5] to obtain 

RHS = T-\a; r) f[ e A °^ r ' )dT ' T(a; r ). (42) 



TO 



Moreover, making use of the defining Eq. (23), we get 

RHS = T-\a- r) Q(a; r, r ) P(a, a ; r„) Q(a; r , r ) = T'\a; r) U(a; r). (43) 

The last line is clearly the expression for W given by Eq. (36). 

Finally, we want to express the Wilson loop W in yet another form which we state 

as: 

Theorem 8 The Wilson loop operator defined in Eq. (26) can be expressed as a 
surface integral of the field strength: 

W = f[ e J «o T-^;r')F m ^;r')T{a';r'Wdr' ^ 

TO 

where Fq\ is the 0-1 component of the non-Abelian field strength. 
To prove this theorem, we note that 

^[T-V,t)4,(<t,t)T(<7,t)] = T-\a,r) [d a A (a,r) + [A (a, r), A^a, r)]] T(a, r). 

(45) 

Moreover, 

^{T-\a, t) (L T T(a, r)) T(a, r)} = T~\a, r)d T Ax(o~, t)T(<t, t). (46) 
It then follows that 

|-{T~V, r)[A (a, r) - L r T(a, r)]}T(a, r)} 
= T-\a, r) [£A (a, r) - r) + [A (a, r),A 1 (a, r)]}T(a, r) 

= T- 1 (a 1 r)F 01 (a,T)T(a,r). (47) 

The last step follows from the definition of the field strength in terms of the connection 
given above 

F 01 ■= foM<r,r) - ^A 1 {a,r) + [A (a, r), A x {a, r)\. (48) 
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Integrating Eq. fl4"T|) with respect to a, we get 

T-\a, r)[A (a, r) - L T T(a, r)]}T(a, r) 
= f T- x {(r'-T')F 1Q {a'-T')T{a'-T')da'dT'. (49) 

We thus arrive at the surface integral representation of the Wilson loop operator: 

W = f[ J°o r-V;OiW;rW;rW<iV (5Q) 



TO 



We note that in this expression the ordering of the operators is defined with respect 
to r whereas a is a parameter. Recalling the antisymmetry of the components of the 
field strength, we can rewrite this expression in terms of path ordered exponentials 
familiar from the physics literature: 

W = Vt6 \ / E da^T-Ha;r)F ab (a;r)T(a;r) ^ ^ 

where da ab is the area element of the 2-surface. Despite appearances, it must be 
remembered that a and r play very different roles in this expression. 

5 A Second Proof 

To illustrate the power and the flexibility of the product integral formalism, we give 
here a variant of the previous proof for the non-Abelian Stokes theorem. This time 
the proof makes essential use of the non-trivial Theorem 6. We start with the form 
of W given in Eq. (36) and take its derivatives with respect to r: 

dW 

—— = d T Q ((T ;t,t )P (a,a ;T)Q(cT]T,T )P(a,a ]T ) + 

+Q~ 1 (a ; t, r )5 r P _1 (cr, a ; r)Q(a; r, r )P(cr, a ; r ) + 
+Q _1 (c r o; t, T )P~ l (cr, a ; r)d T Q(a; r, T )P(a, cr ; r ). (52) 

Here, we have made use of the fact that P(a, a ; r ) is independent of r. As a 
preparation for the use of Theorem |3], we start with Eq. (13) for W and use Theorem 

E 

L T W = ^W- 1 = T- >; t) [A (a; r) - P(a, a ; r)A (a ; r)P-\a, cr ; r) - 

-d T P(a, a ; t)P-\<t, a ; r)] T(a; r). (53) 

Now we can use Theorem 6 to evaluate the derivative of the product integral with 
respect to the parameter r: 

d T P(a,a ;r) = f da'P(a,a' ] T)d T A 1 (a';r)P(a',a ;r). (54) 
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Then, after some simple manipulations using the defining equations for the various 
terms in Eq. (53), we get: 

T-\a; T)d T P(T)p-\T)T(a- t) = r da'T-\a'- T)d T A x {a'\ r)T(a'; r). (55) 

JtT 

Using Theorem || and the fact that P(cr , a ; r) = 1, we can write the rest of Eq. (53) 
as an integral too: 

r~V; r)[A (a; r) - P(a, a ; r)A (a Q - r)p-\a, a Q ; r))T{a; r) = 
= Q- 1 (a ;T,T )[p- 1 (a,(r ;T)A (a;T)P((r,ao;T) - A (a )]Q(a ; r, r ) = 
= /* da' P- V, a ; r){d T A (a', r) + [A (a', r), Al(a', r)])P(a', a ; r). (56) 

Combining Eqs. (53), (55), and (56), we obtain: 

dW 1 r 

L T W = —W~ l = da'T-\a\T)F m {a\T)T{a\r). (57) 

UT J (To 

Using Theorem 3, we are immediately led to Eq. (50) which was obtained by the 
previous method of proof. 

There are two reasons for the relative simplicity of this proof over the one which 
was given in the previous section. One is due to the use of differentiation with respect 
to a parameter according the Theorem 6. The other is due to the use of Eq. (13) and 
Theorem 3 for the composite operator W. In the first proof, the use of this theorem 
for W was not assumed. Its justification for using it in the second proof lies in the 
composition law for product integrals given by Eq. (11). 



6 Convergence Issues 



The definitions of Wilson lines and Wilson loops as currently conceived in the physics 
literature involve exponentials of operators. The standard method of making sense 
out of such exponential operators in the physics literature is through their power 
series expansion: 

A(x)dx ] , (5£ 



Ve J h a A(x)dx = y }_ v ( f b 
n=0 n\ \Ja 



where a typical path ordered term in the sum has the form 

— V{\ A{x)dx) n := / dx 1 dx 2 ... dx n A(x 1 )A(x 2 )...A(x n ). (59) 

TV. Ja Ja Ja Ja 

Without additional specifications, such a power series expansion is purely formal, 
and it is not clear a 'priori that the series fl58|) is well defined and convergent. In- 
deed, in previous attempts [fLOl - f20fl at proving the non-abelian Stokes theorem, the 
convergence of such series has been taken more or less for granted. 
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One important advantage of our product integral approach is that, without the 
need for further input, we can show the convergence of the path ordered exponentials 
precise, without any further effort by using the established properties of product 
integrals. They will enable us to prove that the series of partial sums converges 
uniformly to the product integral. The proof is contained, as a special case, in the 
following two theorems valid for all product integrals. The detailed proofs of these 
theorems are given in reference |§. 

Theorem 9 Given the continuous function A : [a, b] — ► C nxn , and given x, y £ [a, b], 
let L(x,y) = fx \\A(s)\\ds. Also let J (x,y) = I, and for n > 1 define iteratively 
J n (x,y) := J% A(s)J n -i(s,y)ds. Then for any n>0 the following holds: 

y n 1 

\\Ee A ^ ds -J2J k (x,y)\\ < ^^\L(x,y)\ n+1 e L ^ b \ (60) 

This estimate is uniform for all x, y in the interval [a, b]. One of the consequences of 
this estimate is the content of the next theorem. 

Theorem 10 With A and J/~(x,y) as in Theorem 9, we have, in the same notation, 

y oo 

Y[ e ^ = ^J k (x,y). (61) 

The series on the right hand side of this expression converges uniformly for any x, y £ 
[a,b]. 

To give a flavor of the proofs involving product integrals we will include the proof 
of this theorem. According to Theorem |S| we have: 

y n i 

|| JJe*M* - £ J k (x,y)\\ < M r ——\L{b,a)\ n+l — 0, asn - oo (62) 

This follows immediately since, e.g., by Stirling's formula, the asymptotic behavior 
of the factorial function is roughly n\ f=a n n as n — > oo. Therefore, 2_ « (-1 — > as 
n — > oo. q.e.d. 



7 Physical Observables 

From the Wilson loop operator, one can obtain physical quantities in a variety of 
ways. The most familiar one is its trace which gives the c-number Wilson loop or the 
Wilson loop observable: 

trW R (C) = trVe \k<^ h T-\^r)F ab {a-,r)Ti,a-,r) ^ 

^or example, to estimate a general term in the expansion (|58|) to show the norm convergence, we 
don't need to restate the elementary fact that the space of matrices of given rank over the complex 
numbers form a Banach space. 
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The subscript R in this expression refers to the particular representation of the gen- 
erators. 

Another invariant associated with an operator (a matrix) is its determinant. From 
its product integral representation, the determinant of the Wilson loop operator is 
given by 

det W = e tT ^ is ^ b T-Ha; T )F ab (a;r)T(a; T ) ^ 

After some straight forward manipulations, this can be expressed in the form 

detW = e^^ trdaa " Fab{(7 ' T) . (65) 

The generators of simple Lie groups can be represented by traceless matrices so that 
for these groups tr = 0, indicating that det W — 1. This is not surprising since the 
Wilson loop operator is a group element, and for group elements with determinant 
one this result follows trivially. For non-simple groups such as U(l) and the products 
thereof the trace reduces to the trace of commuting elements of the algebra with 
non-zero trace. The corresponding subgroup is commutative, there is no ordering 
problem, and the surface representation of Wilson loop operator reduces to that of 
the (abelian) Stokes Theorem. 



8 The gauge transforms of Wilson lines and Wil- 
son loops 

Under a gauge transformation, the components of the connection, i.e. the gauge 



potentials, transform according to |24 



Ap(x) — > g(x)A^(x)g 1 (x) - g(x)d^g(x) \ (66) 

The components of the field strength (curvature) transform covariantly: 

F^{x) — > g{x)F llv {x)g-\x). (67) 

Using the product integral formalism, we want to derive the effect of these gauge 
transformations on Wilson lines and Wilson loops. 

Let us start with the Wilson line defined by Eq. (22). Under the gauge transfor- 
mation flBTf ) this quantity transforms as 

(7 (7 

ll„i.^\_„l„i. rr \a „-l/V/., 



pfa aQ . r) = TT e M'>'rr)da' > -Q e b(<T , ;r)A 1 (<7';T) 9 - 1 ( ( 7';r)-s(<T';r)9 CT g- 1 (<T';r)]d ( 7'_ 



By Eq. (13), g(a; r)(9 (7 (7 _1 (cr; r) = —L a g(a;T). Thus, we have for the gauge trans- 
formed Wilson line 

j{ e \9(<T';T)A x {a'-,T)g- 1 {a'iT)+L„g{a'-,T)}d*T l _ ^gg^ 
OX) 
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Moreover, we use Theorem 4 and recall from Theorem 3 that Y\° Q e L ^9{^';r)da' _ 
g(a; r)g~ 1 (a ] r). Then, the gauge transform of P(a, a ; r) will take the form 

g(a; r)g-\a Q] t) f[ e^iT)^';^" W). (70) 

Finally, using Theorem 7 and 8 one can readily see that the constant terms in the 
exponents can be factored out from the product integral so that we get 

P(a, a ; r) — ► g(a; r)P(a, a ; T)g~ 1 {a \ r). (71) 

In the physicist's notation, the result can be stated as 



Thus, we have an unambiguous proof of how the Wilson line transforms under gauge 
transformations. This is of course consistent with the role of the Wilson line as a 
parallel transport operator. For a closed path, the points a and b coincide. As a 
result, the corresponding Wilson loop operator transforms gauge covariantly. 

For consistency, we expect that the surface integral representation of the Wilson 
loop also transforms covariantly under gauge transformations. To show this explicitly, 
we note from Eq. (50) that in this case we need to know how the operator T(a, r) 
transforms under gauge transformations. To this end, we note that the Wilson line 
Q(<7] r, T ) given by Eq. (23) transforms as 

Q{°\t,t q ) = j[e*>WW ^(a;r)Q(a;T,To)<rV;To). (73) 

TO 

The transform of the composite Wilson line T(<j,t) given by Eq. (31) follows imme- 
diately: 

T(a; t) = P(a, a ; r) Q(a ; r, r ) — > g(a; r)T(a; r)g-\a o] r ). (74) 

As expected from the composition rule given by Eq. (11), the product of two Wilson 
lines transforms as a Wilson line. 

From the above results, it is straight forward to show that the surface integral 
representation of Wilson loop transforms as 

w ^ e 9(^o;ro)( T-\a'y)F 10 (a'y)T{a'-y)dt^g-\a -, n) ) ^ 

TO 

As in the case of Wilson line, the constant factors in the exponent factorize, so that 
under gauge transformations the Wilson loop transforms covariantly, i.e., 

W — > #(a ;r ) jjy-o g (a ;r ). (76) 

TO 
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This result strengthens our confidence in the self consistency of our formalism. In the 
physics notation, this transformation law takes the form 

where a is a point on the loop C. 

An important consequence of the gauge covariance of the Wilson loop operator 
is that the Wilson loop observable given by Eq. (63) is gauge invariant. Since 
this observable is the trace of the Wilson loop operator, the result follows from the 
invariance of trace of a product of operators under cyclic permutation of the operators. 



9 Concluding Remarks 

The identification of Wilson lines and Wilson loops of non-Abelian gauge theories 
with product integrals allows for the possibility of extracting physical consequences 
from these objects in a consistent and mathematically well defined manner. Although 



many of the properties of Wilson lines and Wilson loops have been discussed [10 



ZD| from various, more intuitive, points of view, there are two issues associated with 



these operators with respect to which their product integral representations have a 
decided advantage. One is the existence issue discussed in Section 6, and the other is 



the supersymmetric generalization of these notions [26]. We are optimistic that the 



present work will help fill the gap in connection with these as well as other issues. 
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number DOE-FGO2-84ER40153. The hospitality of Aspen Center for Physics in the 
Summer of 1998 is also gratefully acknowledged. We are also grateful to Dr. M. 
Awada for valuable input at the initial stages of this work. 



18 



References 

[1] K.G. Wilson, "Confinement Of Quarks," Phys. Rev. D10, 2445 (1974). 

[2] A.M. Polyakov, "Compact Gauge Fields And The Infrared Catastrophe," Phys. 
Lett. 59B, 82 (1975). 

[3] E. Witten, "Topological Quantum Field Theory," Commun. Math. Phys. 117, 
353 (1988). 

[4] M.B. Green, J.H. Schwarz and E. Witten, Superstring Theory, Vol. 1,2. 
Cambridge Monographs On Mathematical Physics. 

J.G. Polchinski, String Theory, vols. I and II, 
Cambridge Monographs on Mathematical Physics, 1998. 

[5] S. Mandelstam, "Gauge Independent Formulation of Electrodynamics," Ann. 
Phys. (N.Y.) 19, 1 (1962). 

[6] T.T. Wu and C.N. Yang, "Concept Of Non-integrable Phase Factors And Global 
Formulation Of Gauge Fields," Phys. Rev. D12, 3845 (1975). 

[7] J. Maldacena, "The Large N Limit of Superconformal Field Theories and Super- 
gravity," Adv. Theor. Math. Phys. 2, 231 (1998), |hep-th/9711200| . 



E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2, 
253 (1998) |hep-th/98021"50| . 



S.S. Gubser, I.R. Klebanov and A.M. Polyakov, "Gauge theory correlators from 
noncritical string theory," Phys. Lett. B428, 105 (1998) |hep-th/ 9802109 . 



E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in 
gauge theories," Adv. Theor. Math. Phys. 2, 505 (1998) |hep-th/9803T3T . 



[8] J. Maldacena, "Wilson Loops in large N Field Theories," Phys. Rev. Lett. 80, 
4859 (1998). [hep-th/ 98030021 

S. Rey and J. Yee, "Macroscopic strings as heavy quarks in large N gauge theory 



and anti-de Sitter superg ravity," |hep-th/980300l 



[9] J.D. Dollard, C.N. Friedman, Product Integration, Addison Wesley, 1979. 

[10] M.B. Halpern, "Field Strength And Dual Variable Formulations Of Gauge The- 
ory," Phys. Rev. D19, 517 (1979). 

[11] Y. Aref'eva, "Non-Abelian Stokes Formula," Theor. Math. Phys. 43, 353 (1980). 

[12] P.M. Fishbane, S. Gasiorowicz and P. Kaus, "Stokes' Theorems For Non-abelian 
Fields," Phys. Rev. D24, 2324 (1981). 



19 



[13] N.E. Bralic, "Exact Computation Of Loop Averages In Two-Dimensional Yang- 
Mills Theory," Phys. Rev. D22, 3090 (1980). 

[14] L. Diosi, "Comments On The Non-abelian Stokes Theorem," Phys. Rev. D27, 
2552 (1983). 

[15] Y.A. Simonov, "Cluster Expansion, Non-abelian Stokes Theorem, and Magnetic 
Monopoles," Sov. J. Nucl. Phys. 50, 134 (1989). 

Y.A. Simonov, "Topological Charges And Convergence Of The Cluster Expan- 
sion," Sov. J. Nucl. Phys. 50, 310 (1989). 

[16] B. Broda in Advanced Electro-magnetism, ed. T. Barrett, D. Grimes, World Sci- 
entific, 1995. 

[17] M. Hirayama and S. Matsubara, "Stokes theorem for loop variables of non- 
Abelian gauge field," Prog. Theor. Phys. 99, 691 (1998) |hep-th/9712T20 . 



M. Hirayama, M. Kanno, M. Ueno and H. Yamakoshi, "Non-Abelian Stokes 
theorem for loop variables associated with nontrivial loops," |hep-th / 9806098| . 

[18] D. Diakonov and V. Petrov, "Non-Abelian Stokes theorem and quark - monopole 
interaction," [hep-th/9606104 . 



[19] F.A. Lunev, "Pure bosonic worldline path integral representation for fermionic 
determinants, non-Abelian Stokes theorem, and quasiclassical approximation in 
QCD," Nucl. Phys. B494, 433 (1997) |hep-th/9609166| . 

[20] E. Alvarez, C. Gomez and T. Ortin, "String representation of Wilson loops," 
|hep-th/ 9806075| . 



[21] A.M. Polyakov, "String theory and quark confinement," Nucl. Phys. Proc. Suppl. 
68, 1 (1997) Eiep-th/97110(J3 

A.M. Polyakov, "The Wall of the cave," |hep-th/ 9809057 . 



[22] L. Schlessinger, Math. Ann. 99 (1927) 413. 

[23] K. Bichteler, Integration Theory, LNM 315, Springer Verlag, 1973. 

[24] Our conventions follow closely M.E. Peskin, D.V. Schroeder, An Introduction to 
Quantum Field Theory, Addison- Wesley, 1995. 

[25] A preliminary version of this work was presented at the International Conference 
on Orbis Scientiae 1998, Fort Lauderdale, Florida, December 17-21, 1998. 

[26] R.L. Karp and F. Mansouri, University of Cincinnati preprint UCTP102.99, and 
in preparation. 



20 



